XY frustrated systems: continuous exponents in discontinuous phase transitions. 
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XY frustrated magnets exhibit an unsual critical behavior: they display scaling laws accompanied 
by nonuniversal critical exponents and a negative anomalous dimension. This suggests that they 
undergo weak first order phase transitions. We show that all perturbative approaches that have been 
used to investigate XY frustrated magnets fail to reproduce these features. Using a nonperturbative 
approach based on the concept of effective average action, we are able to account for this nonuniversal 
scaling and to describe qualitatively and, to some extent, quantitatively the physics of these systems. 
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I. INTRODUCTION 

After twenty-five years of intense activity, the physics 
of XY and Heisenberg frustrated systems is stiU the sub- 
ject of a great controversy concerning, in particular, the 
nature of their phase transitions in three dimensions (see 
for instance Ref. 1 for a review). On the one hand, a 
recent high-order perturbative calculations*^ predicts in 
both cases a stable fixed point in three dimensions and, 
thus, a second order phase transition. On the other hand, 
a nonperturbative approach, the effective average action 
method, based on a Wilson-like Exact Renormalization 
Group (ERG) equation, leads to first order transitionsi. 
Actually, it turns out that, in the Heisenberg case, these 
two theoretical approaches are almost equivalent from 
the experimental viewpoint (see however Ref.l5jj) . Indeed, 
within the ERG approach, the transitions are found to be 
weakly of first order and characterized by very large corre- 
lation lengths and pseudo-scaling associated with pseudo- 
critical exponents close to the exponents obtained within 
the perturbative approach. This occurence of pseudo- 
scaling and quasi-universality has been explained within 
ERG approaches by the presence a local minimum in the 
speed of the flow^, related to the presence of a complex 
fixed point with small imaginary parts, called pseudo- 
fixed point^. 

XY frustrated magnets are rather different from this 
point of view since their nonperturbative RG flows dis- 
play neither a fixed point nor a minimum. We show 
in this article that they nevertheless generically exhibit 
large correlation lengths at the transition and thus, 
pseudo-scaling, but now without quasi-universality. More 
precisely, we show that quantities like correlation length 
and magnetization behave as powers of the reduced tem- 
perature on several decades. A central aspect of our ap- 
proach is that, although the RG flow displays neither a 
fixed point nor a minimum, it remains sufficiently slow 
in a large domain in coupling constant space to produce 



generically large correlation lengths and scaling behav- 
iors. We argue that our approach allows to account for 
the striking properties of the XY frustrated magnets like 
the XY Stacked Triangular Antiferromagnets (STA) such 
as CsMnBrg, CsNiCla, GsMnlg, CsCuGlg, as well as XY 
helimagnets such as Ho, Dy and Tb, which display scal- 
ing at the transition without any evidence of universal- 
ity. Our conclusions are in marked contrast with those 
drawn from the perturbative approach of Pelissetto et al. 
^■^ which leads to predict a second order phase transition 
for XY frustrated magnets. 



II. THE STA MODEL AND ITS LONG 
DISTANCE EFFECTIVE HAMILTONIAN 

The prototype of XY frustrated systems is given by the 
STA model. It consists of spins located on the sites of 
stacked planar triangular lattices. Its hamiltonian reads: 



H 
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where the Si are two-component vectors and the sum 
runs on all pairs of nearest neighbors. The spins inter- 
act antiferromagnetically inside the planes and either fer- 
romagnetically or antiferromagnetically between planes, 
the nature of this last interaction being irrelevant to 
the long distance physics. Due to the intra-plane an- 
tiferromagnetic interactions the system is geometrically 
frustrated and the spins exhibit a 120° structure in the 
ground state (see FIG. ^a). As i/ is invariant under 
rotation, other ground states can be built by rotating 
simultaneously all the spins. 

Let us describe the symmetry breaking scheme of 
the STA model in the continuum limit. In the high- 
temperature phase, the hamiltonian (Q is invariant un- 
der the 5*0(2) X Z2 group acting in the spin space and the 
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FIG. 1: The ground state configurations a) of tiie spins on tlie 
triangular lattice and b) of the order parameter made of two 
orthonormal vectors. The three-dimensional structure of the 
ground state is obtained by piling these planar configurations. 



becomes ferromagnetic. It is thus trivial to derive the 
effective low-energy hamiltonian relevant to the study of 
the critical physics which writes: 



H = -- /"d''xTr(a*$(x).a$(x)) 



(5) 



where *$ denotes the transpose of $. 

It is convenient to consider, in the following, a general- 
ization of the models ^ and © to A^-component spins. 
The order parameter consists in this case in a A^ x 2 ma- 
trix and the symmetry-breaking scheme is thus given by 
0{N)xO{2) -> 0{N~2)xO{2)^^ Frustrated magnets 
thus correspond to a symmetry breaking scheme isomor- 
phic to 0{N) —^ 0{N-2) that radically differs from that 
of the usual vectorial model which is 0{N) — + 0{N — 1). 
The matrix nature of the order parameter together with 
the symmetry breaking scheme led naturally in the 70's 
to the hypothesis of a new universality class i2iS — the 
"chiral" universality class — gathering all materials sup- 
posed to be described by the hamiltonian (^): STA and 
helimagnets. As we now show, examining the current 
state of the experimental and numerical data, there is, 
in fact, no clear indication of universality in the critical 
behavior of XY frustrated magnets. 



III. THE EXPERIMENTAL AND NUMERICAL 
CONTEXT 



0{2) group associated to the symmetries of the triangu- 
lar lattice^. In the low-temperature phase, the residual 
symmetries are given by the group 0(2)^j^^^„ which is a 
combination of the group acting in spin space and of the 
lattice group. The symmetry breaking scheme is given 
byM: 



G = 0(2) X 30(2) X Za ^ i? = 0(2)^^ 



lag 



(2) 



and thus consists in a fully broken SO (2) x Z2 group. The 
Z2 degrees of freedom are known as chirality variables. 

Due to the 120° structure, the local magnetization, 
defined on each elementary plaquette as: 



S = 5i + 52 + 53 



(3) 



vanishes in the ground state and cannot constitute the 
order parameter. In fact, as in the case of colinear antifer- 
romagnets, one has to build the analogue of a staggered 
magnetization. It is given by a pair of two-component 
vectors (j>i and 02 — defined at the center x of each el- 
ementary cell of the triangular lattice — that are or- 
thonormal in the ground stat o'^1^1^ (see FlC^b). They 
can be conveniently gathered into a square matrix: 



$(X) = (01(X),<^2(X)) 



(4) 



Once the model is formulated in terms of the order 
parameter, the interaction, originally antiferromagnetic. 



A. The experimental situation 

Two kinds of materials are supposed to undergo a 
phase transition corresponding to the symmetry breaking 
scheme described above: the STA — CsMnBra, CsNiCla, 
CsMnIa, CsCuClg — (see Ref.Ofor RbMnBrg) and the 
helimagnets: Ho, Dy and Tb. The corresponding critical 
exponents are given in Table Q] 

Note first that, concerning all these data, only one er- 
ror bar is quoted in the literature, which merges sys- 
tematic and statistical errors. We start by making the 
hypothesis that these error bars have a purely statistical 
origin. Under this assumption, we have computed the 
— weighted — average values of the exponents and their 
error bars. This is the meaning of the numbers we give in 
the following. This hypothesis is however too naive, and 
we have checked that, if we attribute a large part of the 
error bars quoted in Table ^ to systematic bias — typi- 
cally 0.1 for /? and 0.2 for 1/ — , our conclusions still hold. 
We also make the standard assumptions that the mea- 
sured exponents govern the leading scaling behavior, i.e. 
the determination of the critical exponents is not signif- 
icantly affected by corrections to scaling. This is gener- 
ically assumed in magnetic materials where corrections 
to scaling are never needed to reproduce the theoretical 
results in the range of reduced temperature reachable in 
experiments . This is different for fluids where the scal- 
ing domain can be very large. Moreover, since the error 



CsMnBra 


Q=0.39(9), 0.40(5), 0.44(5) 

/3=0.21(1), 0.21(2), 0.22(2), 0.24(2), 0.25(1) 

7=1.01(8), 1.10(5); 1^=0.54(3), 0.57(3) 


CsNiCls 


a=0.342(5), 0.37(6), 0.37(8); /3=0.243(5) 


CsMnIs 


a=0.34(6) 


CsCuCla 


a=0.35(5); /?=0.23(2), 0.24(2), 0.25(2) 


Tb 


a=0.20(3); /3=0.21(2), 0.23(4); iy=0.53 


Ho 


/3=0.30(10), 0.37(10), 0.39(3), 

0.39(2), 0.39(4), 0.39(4), 0.41(4) 

7=1.14(10), 1.24(15); i/=0.54(4), 0.57(4) 


Dy 


/3=0.38(2), 0.39(1), 0.39too2 
7=1.05(7); !/=0.57(5) 


STA 
Monte Carlo 


a= 0.34(6), 0.43(10), 0.46(10) 

/3= 0.24(2), 0.253(10); 7=1.03(4), 1.13(5) 

i/=0.48(2), 0.50(1), 0.54(2) 


six-loop 


a=0.29(9); /3=0.31(2); 7=1.10(4); !/=0.57(3) 



TABLE I: Critical exponents of the XY frustrated models, 
from Refs. 2,11-15 and references therein. For CsCuCla the 
transition has been found of first order and the exponents 
mentioned here hold only for a reduced temperature larger 
that 5.10"^ (see Ref. 16). 



bars in frustrated systems are much larger than in the 
usual ferromagnetic systems — by a factor five to ten, 
see Table Q] — neglecting corrections to scaling should 
not bias significantly our analysis. 

Under these assumptions we can analyze the data. We 
find that there are three striking facts: 
i) there are two groups of incompatible exponents. The 
average value of (3, the best measured exponent, for 
CsMnBrg, CsNiCla and Tb — called group 1 — is given 
by /3 ~ 0.23. It is incompatible with that of Ho and Dy 
— group 2 — which is /3 ~ 0.39 (see TableQlfor details). 
Note that for CsCuCla, whose exponents are compatible 
with those of group 1, the transition has been found to 
be very weakly of first order^^. 

ii) the exponents vary much from compound to compound 
in group 1. For instance, the values of a for CsNiCls and 
CsMnBra are only marginally compatible. 
Hi) the anomalous dimension rj is significantly negative 
for group 1. For CsMnBra, the value of r] determined 
by the scaling relation ry = 2(3/v — 1 with /3 = 0.227(6) 
and V = 0.555(21) is ry = -0.182(38). The inclusion of 
the data coming from CsNiCla and Tb does not change 
qualitatively this conclusion. 

Several conclusions follow from the analysis of the 
data. From point i), it appears that materials that are 
supposed to belong to the same universality class dif- 
fer as for their critical behavior. There are essentially 
three ways to explain this. In the first one, the two sets 
of exponents correspond to two true second order phase 
transitions, each one being described by a fixed point. In 
the second, one set corresponds to a true second order 
transition and the other to pseudo-critical exponents as- 
sociated to weakly first order transitions. In the third, 
all transitions are weakly of first order. 

The first scenario can be ruled out since rj is negative 
for group 1 (point Hi)) while it cannot be so in a second 



order phase transition when the underlying field theory 
is a Ginzburg-Landau (^"'-like theorjiH, as it is the case 
herei. The transition undergone by CsMnBra, CsNiCla 
and Tb is therefore very likely not continuous but weakly 
of first order. This would explain the lack of universality 
for the exponents of group 1 (point ii)). 

In the second scenario, the materials of group 2 un- 
dergo a second order phase transition — 77 is found pos- 
itive there — while those of group 1 all undergo weakly 
first order transitions with pseudo-scaling and pseudo- 
critical exponents. Note that although this scenario can- 
not be excluded, it is quite unnatural in terms of the 
usual picture of a second order phase transition. Indeed, 
it would imply a fine-tuning of the microscopic coupling 
constants — i.e. of the initial conditions of the flow — 
for the materials of group 1 in such a way that they lie 
out of, but very close to, the border of the basin of at- 
traction of the fixed point governing the critical behavior 
of materials of group 2. 

The third scenario, that of generic weak first order 
behaviors for the two groups of materials, seems even 
more unnatural, at least in the usual explanation of weak 
first order phase transitions. 

Actually, we shall provide arguments in favor of this 
last scenario. Also, as we shall see in the framework 
of the effective average action, the generic character of 
pseudo-scaling in this scenario has a natural explanation, 
not relying on the concept of fixed point. Then, no fine- 
tuning of parameters is required to explain generically 
weak first order behaviors in frustrated systems. 



B. The numerical situation 

There are no convincing numerical data concerning he- 
limagnets. For the STA system, three different versions 
have been simulated: 

1) the STA itself l^ii^iSa. 

2) The STAR (Ref-IH) — with R for rigid — which 
consists in a STA where the fluctuations of the spins 
around their ground state 120° structure have been 
frozen. This is realized by imposing the rigidity con- 
straint S = S'l -|- S'2 -l- S'3 = at all temperatures. 

3) A discretized version of the hamiltonian Q, called 
the Stiefel V2,2 model'-^. There, one considers a system 
of dihedrals interacting ferromagnetically, which is rep- 
resented on FlC^b. 

At this stage, we emphasize that the rigidity constraint 
E = which is imposed in the STAR, as well as the 
formal manipulations leading to the Stiefel V2,2 rnodel 
affects only the massive — non critical — modes. Thus, 
all the STA, STAR and Stiefel models have the same 
critical modes, the same symmetries and thus the same 
order parameter. One thus could expect a priori that 
they all exhibit the same critical behavior. 

For the STA system, scaling laws are foundiSiiSiSS 
so that a second order behavior could be inferred. 
The STAR and V2.2 models both undergo first order 



transitions^^. Therefore, by changing microscopic details 
to go from the STA model to the STAR or V2,2 models, 
the nature of transitions appears to change drastically. 
This situation indicates that if STA undergoes a genuine 
second order phase transition, the critical behavior of 
frustrated magnets in general is characterized by a low 
degree of universality, a conclusion already drawn from 
the experimental situation. 

With these behaviors one is brought back to the two 
last scenarios proposed in the previous section: i) the be- 
havior of the STA system is controlled by a fixed point 
while the STAR and Stiefel models lie outside its basin of 
attraction ii) all systems undergo first order phase tran- 
sitions. 

In fact, as shown in Ref. ll^l using the two scaling rela- 
tions ?7 = 2(3/ ly — 1 and 77 = 2 — j/v, ij is found to be neg- 
ative in STA systems — although less significantly than 
in experiments — for all simulations where these calcu- 
lations can be performed. One can thus suspect a (weak) 
first order behavior even for the STA system. This hy- 
pothesis is strengthened by a recent work of Itakura who 
has employed Monte Carlo RG techniques in order to in- 
vestigate the critical behavior of both the STA system 
and the Stiefel model—. Using systems with lattice sizes 
up to 126 X 144 X 126, he has provided evidences for weak 
first order behaviors. 

Let us draw first conclusions from the experimental 
and numerical situations. It appears that the critical 
physics of frustrated magnets cannot be explained in 
terms of a single — universal — second order phase tran- 
sition. A careful analysis of the experimental and numer- 
ical data seems to indicate that a whole class of materi- 
als undergo (weak) first order phase transitions. At this 
stage, no conclusion can be drawn about the existence 
or absence of a true fixed point controlling the physics of 
some realizations of frustrated magnets. To clarify this 
issue, we now present the theoretical situation. 



IV. THE THEORETICAL SITUATION 

The early RG studies of the STA and helimagnets — 
and its generalization to iV-component spins — was per- 
formed in a double expansion in coupling constant and in 
e = 4— d on the Ginzburg-Landau- Wilson (GLW) version 
of the model in Refs. 7 9 2 2i23i It has appeared that, for 
a given dimension d, there exists a critical number of spin 
component, called Nc{d), above which the transition is 
of second order and below which it is of first order. Nat- 
urally, a great theoretical challenge in the study of frus- 
trated magnets, has been the determination oi Nc{d). Its 
value has been determined within perturbative computa- 
tion at three- loop order^'*: 

Na{4: - e) = 21.8 - 23.4e -I- 7.096^ -I- 0{e^) . (6) 

Unfortunately, this series is not well behaved since the co- 
efficients are not decreasing fast. It has been conjectured 



by Pelissetto et aliS4 that Nc{2) = 2. Using this con- 
jecture, these authors have reexpressed (O in the form: 



7Vc(4-e) = 2+(2-e)(9.9-6.77e + 0.16e2) + O(e3) . (7) 

The coefficients of this expression are now rapidly de- 
creasing so that it can be used to estimate Nc{d). For 
d = 3 it provides Nc{3) = 5.3 and leads to the conclu- 
sion that the transition is of first order in the relevant 
Heisenberg and XY cases. 

In agreement with this result, the perturbative ap- 
proaches performed at three loops, either in 4 — e or 
directly in three dimensions, lead to a first order phase 
transition for XY systems with a Nc{3) given respectively 
by NciS) = 3.91 (Ref.Hi) and N^S) = 3.39 (Ref.26). 
However, according to the authors, these computations 
are not well converged. It is only recently that a six-loop 
calculation has been performedS*^ directly in three di- 
mensions which is claimed to be converged in the Heisen- 
berg and XY cases. Note that, for values of TV between 
iV ~ 5 and N :^ 7, the resummation procedures do not 
lead to converged results, forbidding the authors to com- 
pute Nci'i) in this way. For iV = 2 and iV = 3 a fixed 
point is found. The exponents associated to the N = 2 
case are given in Table Q] Note that 7 and ly compare 
reasonably well with the experimental data of group 1. 
However, as we already stressed, the existence of a fixed 
point implies 77 > — rj — 0.08 in Ref.^ — and is thus 
incompatible with the negative value of 77 found for the 
group 1. Moreover, the value p = 0.31(2) found in Ref.jj 
is far — 4 standard deviations — from the average exper- 
imental value (3 — 0.237(6) for group 1 and also far — 3.7 
standard deviations — from that obtained from group 2, 
Ho and Dy: f3 — 0.388(7). It is thus incompatible with 
the two sets of experimental values. This point strongly 
suggests that the six-loop fixed point neither describes 
the physics of materials belonging to group 2 that, in 
the simplest hypothesis, should also undergo a first order 
phase transition. 

The preceding discussion does not rule out the exis- 
tence of the fixed point found in Ref.S. This just shows 
that, if it exists, it must have a very small basin of attrac- 
tion, and that the initial conditions corresponding to the 
STA and helimagnets lie out of it. In fact, as we argue in 
the following, this fixed point probably does not exist at 
all so that we expect that all transitions are of (possibly 
very weak) first order. 



V. THE EXACT RG APPROACHES 

There exists an alternative theoretical approach to the 
perturbative RG calculations which explains well, qual- 
itatively and to some extent quantitatively, all the pre- 
ceding facts. It relies on the Wilsonian RG approach to 
critical phenomena, based on the concept of block spins 
and scale dependent effective theorie322iS&. Although it 
has been originally formulated in terms of hamiltonians, 



its most recent and successful implementation involves 
the effective (average) actionS2i22i2i. In the same way as 
in the original Wilsonian approach, one constructs an ef- 
fective action, noted Fj,, that only includes high-energy 
fluctuations — with momenta q^ > k^ — of the micro- 
scopic system. At the lattice scale fc = A = a^^, Tk 
corresponds to the classical Hamiltonian H since no fluc- 
tuation has been taken into account. When the running 
scale k is lowered, Tk includes more and more low-energy 
fluctuations. Finally, when the running scale is lowered 
to k — 0, all fluctuations have been integrated out and 
one recovers the usual effective action or Gibbs free en- 
ergy F. To summarize one has; 



(8) 



Note also that the original hamiltonian depends on the 
original spins while the effective action — at fc = — 
is a function of the order parameter. At an intermediate 
scale k, Tk is a function of an average order parameter at 
scale k noted (jikil) — or more simply, 4>{q) — that only 
includes the fluctuations with momenta q^ > k"^. Thus 
Tk has the meaning of a free energy at scale k. 

At a generic intermediate scale fc, F^ is given as the 
solution of an exact equation that governs its evolution 
with the running scaled: 




dTk[4>] 



dt =1^' 



(ri^^M 



Rk 



-1 dRk 



dt 



(9) 



where t — In (fc/A) and the trace has to be understood as 
a momentum integral as well as a summation over inter- 
nal indices. In Eq.®, F^ '[(/)] is the exact field- dependent 
inverse propagator i.e. the second functional derivative 
of Tk . The quantity Rk is the infrared cut-off which sup- 
presses the propagation of modes with momenta q^ < k^. 
A convenient cut-off, that realizes the constraints (jHl, is 
provided byS: 



Rkiq^)^Z{k^-q^)eie-q') 



(10) 



where Z is the fc-dependent field renormalization. 

Ideally, in order to relate the thermodynamical quan- 
titites to the microscopic ones, one should integrate the 
flow equation starting from k — A with the microscopic 
hamiltonian H as an initial condition and decrease fc 
down to zero. However, Eq.© is too complicated to be 
solved exactly and one has to perform approximations. 
To render Eq.© manageable, one truncates the effective 
action Tk [4>\ to deal with a finite number of coupling con- 
stants. The most natural truncation, well suited to the 
study of the long distance physics of a field theory, is to 
perform a derivative expansion"^^ of F^ [</)] . This consists 
in writing an ansatz for Tk as a series in powers of d4>. 
The physical motivation for such an expansion is that 
since the anomalous dimension 77 is small, terms with 
high numbers of derivatives should not drastically affect 
the physics. 



Actually, another truncation is performed. It consists 
in expanding the potential, which involves all powers of 
the 0{N) X 0(2) invariants built out of 0i and 02, in 
powers of the fields. This kind of approximation allows 
to transform the functional equation ^ into a set of or- 
dinary coupled differential equations for the coefficients 
of the expansion. It has been shown during the last ten 
years that low order approximations in the field expan- 
sion give very good results (see Ref.|3J for a review and 
Ref.|3^ for an exhaustive bibliography). The simplest 
such truncation is*: 



Tk = 



d'^xl -Tr(a'$.a$) + -(0i. 

z 4 



4 V2 



(11) 



Let us first discuss the different quantities involved in 
this expression. One recalls that $ is the N x 2 ma- 
trix gathering the A^-component vectors 0i and (f>2 (see 
Eq.Ql). There are two independent 0{N) x 0(2) in- 
variants given hy p ~ Tr*$$ and r ~ iTr(*<I>$)^ - 
i(Tr *$$)2. The set {k,X, ii,Z,uj} denotes the scale- 
dependent coupling constants which parametrize the 
model at this order of the truncation. The first quan- 
tity in Ea. ((TT|l corresponds to the standard kinetic term 
while the third and fourth correspond to the potential 
part. Actually, apart from the second term — called 
the current term — , Ffc in Ea. (|ll|l looks very much like 
the usual Landau-Ginzburg- Wilson action used to study 
perturbatively the critical physics of the 0{N) x 0(2) 
model, up to trivial reparametrizations. There is however 
a fundamental difference since we do not use F^; within 
a weak-coupling perturbative approach. This allows the 
presence of the current term which corresponds to a non- 
standard kinetic term. This term is irrelevant by power 
counting around four dimensions since it is quartic in the 
fields and quadratic in derivatives. However its presence 
is necessary around two dimensions to recover the results 
of the low-temperature approach of the nonlinear sigma 
(NLcr) model since it contributes to the field renormal- 
ization of the Goldstone modes. Being not constrained 
by the usual power counting we include this term in our 
ansatz. Note also that we have considered much richer 
truncations than that given by Eq. Hll|l by putting all the 
terms up to ^^'^ and by adding all terms with four fields 
and two derivatives. This has allowed us to check the 
stability of our results with respect to the field expan- 
sion. 

We do not provide the details of the computation. The 
general technique is given in several publications and its 
implementation on the specific 0{N) x 0(2) model will 
be given in a forthcoming article^. The /? functions for 
the different coupling constants entering in l|ll|l are given 
by: 



cLk, 
It 



dX 



{d-2 + ri)K + Avd 



iz^i(0, 0, Ku) + {N- 2)Zfo(0, 0, 0) + hfoiKX, 0, 0) + (l + 2^) if^infi, 0, 0)- 



+ -;2+>,0,Kc.) 



— =id~4 + 2rj)X + Vd 



2X'^{N - 2)Z^o(0, 0, 0) + A2;^2(0> 0) «:^) + SA^'^oM' 0' 0) + 2(A + 2fj.)^l^Q{Kfi, 0, 0) + 



2i4+df 



+ 4Aw/^+''(0, 0, Kuj) + 4w^/^+''(0, 0, Kw) 



dt 



2- '^ ^ ^ 3(2A + /.)„ , , ^^, 8A + /i „ 



-Zgi(0,0,M , 

K k(/i — A) 



7fo(A.A,0,0) + 



k(A — /i) 



Zfo(«M,0,0)- 



+ /^Zfil^A^, 0, tiuo) + ti{N - 2)4(0, 0, 0) 

^ = 2^ (4 - d)Kw?[(i(0, 0, Ku) + 2k'^lo'^11+'^{Q, 0, kuj) + 2m;52(0, 0, kw) - 4m^i(0, 0, kw)+ 

Gt QjKj 

+ 2(-2 + d)Kujl%{0, 0, 0) + 2m^o(0: 0, kuj) + 2k2a2to^ al'*^' 0, 0) + 4kV^™2,2('«M, 0, 0) + 



(12a) 



(12b) 



(12c) 



(12d) 



+ 4Ka;no2(Oi 0, kuj) - 8KUjnf^{0, 0, kuj) + 4Kwn2o(0, 0, kw) 



KUJ< ii— ^;^^(0, 0, KW) + ^^-— ^/^i(kA, 0, KCJ) + Kw/^+'^(0, 0, KUJ)- 



duj , , „ „ , 4:Vd 

-^d-2 + 2,). + -^.„.^ 2 ^^^^---'' 2 

- 3Kw/2+'i(KA, 0, Kw) + {d- 2)lfo{0, 0, 0) - (d - S)lfo{KX, 0, 0) + 8kA/^i(kA, 0, kuj) + 2kujII+'^{k^i, 0, 0) (12e) 
+ 2kuj{N - 2)Z2+'*(0, 0, 0) I + m'l^{0, 0, Ktj) - m'^^{KX, 0, kw) - 2mfi(0, 0, kuj) + 2to'j'i(kA, 0, kw) + 

+ m2o(0, 0, Ktj) - m2o(KA, 0, kcj) + K^A^m22(KA, 0, 0) + 2k'^ ^'^ 11122 {k^, 0, 0) + 2^^77,02(0, 0, kuj)- 

— 4KujnQ2{KX, 0, Kw) — 4:KUjnii{0, 0, koj) + 8KUjnii{KX, 0, kw) + 2Ka;n2o(05 0, kuj) — 4Kujn2Q{KX, 0, kuj) 



In these equations appear the dimensionless functions loop2iS242^: 
^ni n2j™ni,n2:"ni.n2; cahed threshold functioHS, since 
they govern the decouphng of the massive modes enter- 
ing in the action l|ll|) . They encode the nonperturbative 
content of the flow equations p2|l . They are comphcated 
integrals over momenta and are explicitly given in the 
Appendix. 



VI. RESULTS AND PHYSICAL DISCUSSION 



Checks of the method 



§ - (^ - 4)A + Y^ (4Aa^ + V + \'iN + 4)) 
f^(.-4),+ ^(6A, + A.,^). 



(13) 



This can be easily verified considering the asymptotic 
expressions of the threshold functions that are given in 
the Appendix. 

2) Also, around d — 2, performing a low-temperature 
expansion — corresponding to a large k expansion — of 
our equations and making the change of variables: 



rji ~ 2nK 

ri2 ~ 47rK(l -|- kuj) 



(14) 



We have first proceeded to all possible checks of our 
method by comparing our results with all available data l^^ nonfinear a" m^odefat one Toop 
obtained within the different pertubative approaches. 
Our method fulfills all these checks. 



one recovers the /3-functions found in the framework of 

36. 



1) Around d = 4, we have checked that, in the limit 
of small coupling constant, our equations degenerate in 
those obtained from the weak coupling expansion at one 



/3i = -(d - 2)771 + iV - 2 - ^ 

2rii 



/32 = -(d- 2)772 



N-2 frn 



(15) 



This shows that our method allows to recover the per- 
turbative results both near d — 2 and d = 4. This is not a 
surprise since, as well known, Eq.© has a one- loop struc- 
ture. However, for frustrated magnets, and contrary to 
the 0{N) case, the matching with the results of the NLct 
model is not trivial since it requires to incorporate the 
non trivial current term which is irrelevant around d = 4, 
and thus absent in a LGW approach. 

3) Our method also matches with the leading order 
results of the 1/N expansion of ly and rj. In the 0{N) x 
0(2) model these exponents have also been computed at 
order 1/N'^ in d = 3. They are given bjiSi: 



i/= 1- 



V 



56 



16 1 
4 1 64 1 



0{1/N^) . 



(16) 



We have computed v and 77 for a large range of values 
of N to compare our results with those obtained within 
the 1/N expansion. We meet an excellent agreement — 
better than 1% — for v already for A^ > 10 where the 
1/A^ expansion is reliable. 

4) For A^ = 6, a Monte Carlo simulation has been per- 
formed on the STA model^i. A second order phase tran- 
sition has been found without ambiguity with exponents 
given in Table ^ Our results compare very well with 
the Monte Carlo data. Although the case iV = 6 does 
not correspond to any physical system this constitutes a 
success of our approach from a methodological point of 
view. Let us recall that, in this case, the six loop weak 
coupling calculation is not converged. 





a 


P 


7 


V 


■n 


MC 


-0.100(33) 


0.359(14) 


1.383(36) 


0.700(11) 


0.025(20) 


P.W. 


-0.121 


0.372 


1.377 


0.707 


0.053 



TABLE IL Critical exponents for the iV = 6 STA system. The 
first row corresponds to Monte Carlo data;— and the second 
(P.W.) to the present work. 

5) In agreement with the weak coupling perturbative 
results, we have found by varying A^ in a given dimension 
d that, below a critical value Nc{d) of N, no fixed point 
exist3^(see FIC.EJ. Our value of Nc{d) agrees within ten 
percents for all dimensions with the values obtained from 
Eq.O. In d = 3 we have found A^c(3) = 5.1 which almost 
coincides with that obtained using Eq.ljTI) that leads to 
N,{3) = 5.3. 

All these checks show the consistency between our 
computation and most of the previous theoretical ap- 
proaches. It is in contradiction with the six loop calcu- 
lation of Pelissetto et al.. Note that the result A''c(3) > 3 
docs not exclude a priori the existence of a fixed point 
not analytically related — in A^ and e = A — d — to 
those found above Nc{d). This is, in fact, the position 
advocated by Pelissetto et alM. However, we have nu- 
merically searched for this fixed point with our equations 
without success. This implies three possibilities. The 



Nc 




FIG. 2: The full line represents the curve Nc{d) obtained by 
the three- loop result improved by the constraint Nc{2) — 2, 
Eq.(|7|l. The crosses represent our calculation. 



first one is that our method is not able — in principle 
— to find this kind of fixed point. Let us emphasize 
that, although it is not possible to exclude this case, it 
is improbable that a method which recovers all previous 
results in a convincing way, misses such a fixed point. 
Another possibility is that it has a so small basin of at- 
traction that we have systematically missed it in our nu- 
merical investigation of our equations. This would mean 
that it probably does not play any role in the physics of 
frustrated magnets. The transitions should therefore be 
of first order for almost all real systems. It would then 
be very difficult to explain the occurence of generic weak 
first order phase transitions. Finally, there remains the 
possiblity that the fixed point does not exist at all which 
now appears as the most probable solution. 



B. Scaling without fixed and pseudo- fixed point 

The hypothesis of the absence of fixed point imme- 
diately raises the theoretical challenge to explain the 
occurence of scaling in absence of a fixed point. For 
Heisenberg spins, this question has been addressed by 
Zumbach'^^ using a local potential approximation (LPA) 
of the Polchinski equation'^^ and by the present authors 
beyond LPA in the framework of the effective average 
action method^. It has been realized that, when A^ is 
lowered from A^ > A^c(3) to A^ < A^c(3), ahhough the 
stable fixed point disappears, there is no major change 
in the RG flow. 

This can be understood by considering a domain £ of 
initial conditions of the flow corresponding to all systems 
we are interested in — STA, STAR, Vn^2i real materials, 
etc — and studying the RG trajectories starting in £. 
One finds that there exists a domain V C £ such that 
all trajectories emerging from P are attracted towards a 
small domain TZ in which the flow is very slow. Since the 
flow is very slow in TZ the RG time spent in this region is 
long and, thus, the correlation lengths of the systems in 
V are very large. One therefore partly recovers scaling 
for systems in V (that aborts only for very small reduced 



temperatures). Moreover, the smallness oiTZ ensures the 
existence of (pseudo-)universahty. Consequently TZ mim- 
ics a true fixed point. 

This idea has been formahzed through the concept 
of pseudo-fixed point, corresponding to the point in 
TZ where the flow is the slowest, the minimum of the 
flow^. At this point it has been possible to compute 
(pseudo-)critical exponents characterizing the pseudo- 
scaling (and pseudo-universality) encountered in Heisen- 
berg frustrated spin systems^'^^. 

Within our approach, we confirm the existence, for val- 
ues of TV just below Nc{3), of a minimum of the flow lead- 
ing indeed to pseudo-scaling and quasi-universality (see 
Ref.|j and, for details, Ref._5). However, when N is low- 
ered, the minimum of the flow is less and less pronounced 
and, for some value of N between 2 and 3, it completely 
disappears. Since (pseudo-)scaling is observed in exper- 
iments and numerical simulations in XY systems this 
means that the minimum of the flow does not consti- 
tute the ultimate explanation of scaling in absence of a 
fixed point. At this stage, one reaches the limits of the 
notion of minimum of the flow as the quantity playing 
the role of (pseudo-)fixed point. First, it darkens the im- 
portant fact that the notion relevant to scaling is not the 
existence of a minimum of the flow but that of a whole 
region TZ in which the flow is slow, i.e. the /3 functions 
are small. Put it differently, the existence of a minimum 
of the flow does not guarantee that the flow is slow, i.e. 
that the correlation length is large compared with the 
lattice spacing. Reciprocally, it can happen that the RG 
flow is slow, the correlation length being large, and that 
scaling occurs even in absence of a minimum. Second, 
reducing 7?. to a point, one rules out the possibility to 
test the violation of universality. Clearly, the degree of 
universality is related to the size of TZ. The smaller TZ 
the more universal is the behavior. 

Qualitatively — thanks to continuity arguments — one 
expects that, for N close to Nc{3), all systems in £ exhibit 
pseudo-scaling — T> — £ — and TZ is almost pointlike so 
that the transitions are extremely weakly of first order 
and universality almost holds. When N is sufficiently 
decreased, two phenomena occur. First, T> gets smaller 
than £ and thus the transitions become of strong first or- 
der for all systems that belong to £ but not to V. Second, 
TZ gets wider and thus a whole spectrum of exponents 
is observed and pseudo-universality breaks down. Note 
that these two phenomena are not obligatorily simulta- 
neous so that, for intermediate values of N, scaling still 
holds while pseudo- universality is already significantly vi- 
olated. These two phenomena are observed numerically 
in the Heisenberg and XY systems: STA, STAR and 
Vm.2 models. For TV = 3, all these models display scal- 
ing but with (3 exponents that are almost incompatible 
[/3sTA=q^85(ll)^/3sTAii=0.221(9) and /3y3_, =0.193(4) 



order. 

From a theoretical point of view, the Heisenberg case 
has been already treated in Ref.jj. We now concentrate 
on the XY case. 



C. Our results 

In practice, we numerically integrate the flow equa- 
tions H12I) and compute physical quantities like correla- 
tion length and magnetization as a function of the re- 
duced temperature tr ^ {T — Tc)/Tc. Since one expects 
the behavior of frustrated magnets to be nonuniversal, 
one should study each system independently of the oth- 
ers. Thus, ideally, one should consider as initial condi- 
tions of the RG flow all the microscopic parameters char- 
acterizing a specific lattice system. This program would 
require to identify and to deal with an infinite number of 
coupling constants. This remains a theoretical challenge. 
Rather than doing this, wc have chosen to address the 
question of scaling in absence of a fixed point indepen- 
dently of a given microscopic system, using for F^. a finite 
ansatz similar to, but richer than, Ea. (|ll|l . 

Since our truncations forbid us to relate precisely 
the thermodynamical quantities to the microscopic cou- 
plings, we have chosen, as initial conditions, the simplest 
temperature dependence of the parameters at the scale 
k — K. This consists in fixing all coupling constants to 
temperature-independent values and in taking the usual 
ansatz for the quadratic term of Ffe : 



Kk=A 



bT 



(17) 



(see Ref. klf) ] which means that one probably starts leav- 
ing the pseudo-universal regime. For N = 2, scaling is 
only observed for STA — (3 = 0.24(2) here — while the 
transitions for STAR and V2 2 are found to be of first 



where a and b are parameters that we have varied to test 
the robustness of our conclusions. For each temperature 
we have integrated the flow equations and deduced the t^- 
dependence of the physical quantities around the critical 
temperature. 

Let us now review the main results obtained by the 
integration of the RG flow. 

1) The integration of the flow equations leads generi- 
cally to good power laws in reduced temperature for the 
magnetization, the correlation length (see FIG. OJ and 
the susceptibility. We are thus able to extract pseudo- 
critical exponents varying typically, for /? between 0.25 
and 0.38 and for i' between 0.47 and 0.58. This phe- 
nonomenon holds for a wide domain T> in coupling con- 
stant space. 

2) We easily find initial conditions leading to pseudo- 
exponents close to those of group 2: (3 = 0.38, v = 0.58 
and 7 = 1.13 (see Ho and Dy in Tabled for comparison). 
We have checked that the previous result is quite stable 
to both variations of the microscopic parameters and to 
a change of the T-dependence of the microscopic cou- 
pling constants. This is in agreement with the stability 
of P in group 2. In the region of parameters leading to 
this behavior, correlation lengths as large as 5000 lattice 
spacings are found. 




logiol^rl 



FIG. 3: Log-log plots of the magnetization m and of the corre- 
lation length ^ for A'^ = 2 as functions of the reduced temper- 
ature tr for parameters corresponding to materials of group 
1. The straight lines correspond to the best power law fit of 
the data. The power-law behavior observed breaks down for 
small tr- 



3) We also easily find initial conditions leading to 
/3 ~ 0.25, corresponding to group 1. The power laws 
then hold on smaller range of temperature and the criti- 
cal exponents are more sensitive to the determination of 
Tc and to the initial conditions, in agreement with point 
ii) of Section III-A. For such values of /?, we find a range 
of values oi v — 0.47 < v < 0.49 — which is somewhat 
below the value found for CsMnBrs (see Tabled). Also, 
the corresponding rj deduced from the scaling relation 
77 = 2f3/iy — 1 is always positive and, at best, zero. Fi- 
nally, it is very interesting to notice that when we find /3 
of order 0.25 (group 1) we also find correlation lengths at 
the transition of the order of a few hundreds lattice spac- 
ings which coincide rather well with the first size where a 
direct evidence of a first order transition has been seen in 
Monte Carlo RG simulations (lattice sizes around 100)^"'^. 

4) For f3 ~ 0.3 we find critical exponents in good agree- 
ment with those obtained by the six-loop calculation of 
Ref.|3 (see Table Hj). For instance, for /3 = 0.33 we find 
i/ ~ 0.56 and 7 ~ 1.07. 



of XY frustrated systems in three dimensions could un- 
dergo a second order phase transition characterized by 
critical exponents associated with a new universality 
class. We have given convincing arguments that rather 
favor the occurence of — generically weak — first order 
transitions for all XY frustrated magnets with a spread- 
ing of (pseudo-)critical exponents. This is supported by 
experimental and numerical results that do not agree 
with a second order behavior. Moreover we have shown, 
using a nonperturbative approach, that this generic but 
nonuniversal scaling finds a natural explanation in terms 
of slowness and "geometry" of the flow. Our approach 
appears to explain the main puzzling features of XY frus- 
trated magnets. 

We now propose several tests to confirm our approach. 
On the experimental side, more accurate determinations 
of critical exponents could lead to a definitive answer on 
the nature of the transition, at least if no drastically new 
physics emerges (as it could be the case for Helimagnets) . 
In particular, it is important to check that 77 is negative 
for CsMnBr3, and therefore to refine the determination 
of v. It would also be of utmost interest to have more 
precise determinations of a in CsMnBra, CsNiCla and 
CsMnIa which are, up to now, only marginally compat- 
ible. In case they are different, this would corroborate 
the lack of universality that we predict. The existence 
of a continuous spectrum of critical exponents could be 
directly tested by simulating, for example, a family of 
models, extrapolating continuously from STA to STAR. 

On the theoretical side, it would be of interest to 
push the derivative expansion to refine the value of r/ for 
group 1, which is, as usual, overestimated^. This would 
allow to reproduce its observed negative value. More- 
over, it would be interesting to clarify the discrepancy 
between the nonperturbative approach and the six loop 
result. The ability of our approach to reproduce exactly 
the whole set of exponents found by Pelissetto et al. sug- 
gests that their fixed point does not correspond to a true 
fixed point but, in fact, to a region where the flow is very 
slow. One can thus question the convergence of the per- 
turbative result which is not Borel summable. A detailed 
analysis of this problem of convergence could reveal that 
the real fixed point found in the perturbative approach 
is, actually, a complex one. 

Finally, the major characteristics of XY-frustrated 
magnets i.e. the existence of scaling laws with contin- 
uously varying exponents are probably encountered in 
other physical contexts, generically systems with a crit- 
ical value Nc of the number of components of the order 
parameter, separating a true second order behavior and 
a naively first order one. 



VII. CONCLUSION 

On the basis of their specific symmetry breaking 
scheme it has been proposec&^iS that the critical physics 
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equations (|12|l . 



APPENDIX A: THRESHOLD FUNCTIONS 



We discuss in this appendix the different threshold 
functions ^^^.n^, 'm-ti,n2 ^^'^ "'ni,n2 appearing in the flow 



The threshold functions are defined as: 






1 f°° ~ ( 1 



1 poc 



2 7o '\{Pi{y) + wir^{P2{y) + w2Y^ 

- rdyy'^/^^%1 ^^^^^(^) 

2 7o *l(Pi(y)+z«i)"^(P2(2/) + ^2)"^ 



(Ala) 
(Alb) 
(Ale) 



r 



with: 



Piiy) = y{l + r{y) + a) 
P2iy) ^ y{l + r{y)) . 



(A2) 
(A3) 



In all the previous expressions y is a dimensionless quan- 
tity: y — q^ /k'^ where q is the momentum variable over 
which the integral in Eq.@ is performed. As for r{y), it 
corresponds to the dimensionless renormalized infrared 
cut-off: 



r{y) 



Rk{q^) _ Rkiyk^) 



Zq^ 



Zk?y 



(A4) 



In Eas. lAlaVlAlcl dt means that only the t-dependence 
of the function R^ is to be considered and not that of the 



J 



coupling constants. Therefore one has: 



dtPM = ^-i^4rPM 



dt dRk 

-yiwiy) + "^yr'iy))- 



(A5) 
(A6) 



Now the threshold functions can be expressed as explicit 
integrals if we compute the operation of dt- To this 
end, it is interesting to notice the equality: dydtPi{y) = 
dtdyPi{y), so that: 

dtdyT = -T]{r + yr') - 2y(2r' + yr") (A7) 

We then get: 






1 [°° 
"ni,«2("'i'^2,a) = ^9 / dy 



,,d/2_ 



(Pi(2/)+wi)«i(F2(y)+W2)"^ 



y (1 + a + 7' + yr'){rjr + 2yr') x 



ni 



n2 



Pi{y)+wi P2iy)+W2 



rj{r + yr') — 2y{2r' + yr") 



™«i.n2(^l'^2,a) 



dyy' 



n-i 



d/2_ 



1 -I- a -I- r -f yr' 



{Piiy)+wi)"^{P2{y)+W2)''^ 



y(l + a + r + yr') (r/r + 2yr') x 



n2 



Piiy) + wi P2{y)+W2 



277 (r + yr') - Ay{2r' + yr") 



(A8) 



(A9) 



(AlO) 
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